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VI. 1 PERFORMANCE 


Aircraft performance modeling requires consideration of propulsion, 
aerodynamics and weight characteristics, and the interaction between 
the subroutines may be quite complex. This is particularly true when a 
large number of realistic physical constraints are included in the per- 
formance specification. Table VI. 1.1 shows the interaction between the 
subroutines as discussed in this volume of the study. 


Table VI. 1.1 



PERFORMANCE SUBROUTINES 

PERFRM (70 cards) 

DEFINES MISSION 

ACCEL ( 70) 

RECTILINEAR SPEED-UP 

CLIMB (250) 

CONSTRAINED CLIMBING FLIGHT 

DLAND (170) 

LANDING MANEUVER 

TAKOFF (350) 

TAKEOFF MANEUVER TO FLAPS-UP 

DERIV ( 90) 

DETERMINES DERVIATIVES AND CHECKS CONSTRAINTS FOR 


TAKOFF ROUTINE 

TAXI ( 25) 

FUEL CONSUMPTION DURING GROUND IDLE 

TURN ( 40) 

CONSTANT ALTITUDE TURN MANEUVER 

X RANGE (340) 

CRUISE RANGE COMPUTATION 

ASPEED ( 90) 

DETERMINING SPEED AT SPECIFIED POWER FOR EQUILIBRIUM 


HORIZONTAL FLIGHT 

RGBAL (120) 

CONTROL ROUTINE WHEN A RANGE OR ENDURANCE IS 


SPECIFIED IN INPUT 


VI. 1.1 Subroutine PERFRM 


This subroutine is called only by MAIN, and it acts principally 
to organize and unify the computations carried out by other performance 
subroutines. The call has the form CALL PERFRM I cqnd' I FLY^ ’ T ^ e 
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mission segment number is I___, which is associated with a specific 

SEG 

subroutine; i.e., an shown in Figure VI. 1.1, the mission segment varies 
as follows: 

Call TAXI 
Call TAKOFF 
Call ACCEL 
Call CLIMB 

Write Performance Data and Advance I 

SEG 

Call XRANGE 
Call DLAND 
Call TURN 
Return 



The other call parameters are 


X COND 


jO, Write data heading and set 1^^ • 1 
1 1 » Initialize weight, altitude, etc., and continue 


! 1, Call all performance subroutines, unless I QC ■ 99 
2, Call TAXI and TAKOFF only 
3, Call TAXI, TAKOFF and CLIMB only 


The other significant input parameters to the subroutine will be 

described as they are needed in the discussion to follow. 

The program logic beqins before statement 5, where I___ is aug- 

SEG 

mented for later use; 


I 


SEG 


*SEG * 1 


(VI. 1.1) 


and the gross weight W nd altitude H are specified in terms of the 

input parameters WT misn , and The various mission segments 

are then analyzed in the order implied by I . The basic order is 

SEG 

taxi, takeoff, climb, cruise, reserve, and landing with acceleration 
between takeoff and climb and climb and cruise. 
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Return 


Figure VI. 1.1 Subroutine PERFRM Flow Chart 
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For I seG ■ 1, the subroutine TAXI requires as input the runway 
altitude H, the initial time S^, and the taxi time, DELTT. It returns 
the parameters 

W - grojs weight after taxi, lb 
WF - fuel burned in taxi, lb 
ST ■ time at end of taxi, sec 

and then if I ? 99, I is augmented. 

DC SEO 

For I * 2, the subroutine TAKOFF is called, which computes the 
SEG 

time history of the take-off maneuver, from ground roll to a specified 
height. The standard maneuver is divided into sub-segments, or "events," 
as follows: 


(1) 

Begin ground roll 




(2) 

Rotate aircraft 




(3) 

Lift off 




(4) 

Begin gear retraction 




(5) 

Distance to an altitude 

of 

35 

ft 

(6) 

Distance to an altitude 

of 

50 

ft 

(?) 

Begin flap retraction 




(8) 

End take-off 





In addition, for multi-engine aircraft, the optional engine out 
takeoff performance can be computed. This consists of the continued 
takeoff performance to an altitude of 50 ft with an engine failed and 
the acceierate-stop distance divided into the events, 

(1) Begin ground roll 

(2) Engine failure 

(3) Remove power 

(4) Apply brakes 

(5) Stop aircraft 

The subroutine ACCEL (SPEED) is called, between mission segments 
to accelerate the aircraft to the speed of the next segment. This subroutine 
models the rectilinear speed-up maneuver at constant altitude, during 
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which weight and speed change monotonicaily, until the velocity reaches 
the value SPEED. SPEED may be specified in knots or Mach number. 

When I„_ is 4, CLIMB is called, with the altitude argument 

5 EG 


h endcl 


cr alt* EH rURN " 0 

h turn* “Surn * 0 


(VI. 1.2) 


EM^g^ is the Mach number during the turn. If turning performance is 

desired it will be computed after the climb otherwise the climb is to 

the cruise altitude CR. . . If no value is input for CR it will be 

ALT ALT 

defaulted to EM CR(J . 

The cruise portion of the mission (I_„_ ■ 6) is analyzed by a call 

SEG 

to subroutine X RANGE ( W FCRU « w rRES ) ' where 


W * fuel used during cruise, lb 

FCRU 

W » reserve fuel, lb 

FRES 


are calculated by X RANGE in terns of three different payloads and three 
different Mach numbers. 

Subsequently, I __ is advanced to 7, for which subroutine DLAND is 

SEG 

called, using the flag I * 99. This subroutine computes the landing 

IaER 

field performance. 

:r the optional turning performance is desired, subroutine TURN is 
called, after the climb segment. This brief subroutine calculates 
aircraft turning performance at constant altitude, Mach number, and load 
factor, and it calculates various characteristics of this steady state 
maneuver. 


VI . i . 2 Subroutine ACCEL 

This subroutine models the thrust and drag characteristics of an 
aircraft while its forward velocity is changed from an initial value 
to a required final value. 

The subroutine begins by calling TPALT, which returns various 
atmospheric properties existing at the initial altitude. The static 
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pressure P q is converted to lb per sq ft, and the speed of sound is 
given in ft/sec as 

S - 49.1 SQRT (T ) (VI. 1.3) 

A O 

The air density is a function of pressure, temperature and gravity, 

RHO - P / [53.32 T G) (VI. 1.4) 

o o 

and these units are in slugs per cu ft. The single calling argument 
of the subroutine is a command Mach number or velocity, SPEED, which 
must be greater than the current Mach number, EM, in order for the 
positive acceleration modeling to be effective. If the calling argument, 
SPEED, is greater than 5.0 the argument is considered to be the command 
velocity in knots. The initial velocity in ft per sec is given by 

V ■ EM * S (VI. 1.5) 

A 

the time-step interval is 0.5 sec, and the velocity difference to be 

made up is 


DE1 V " X MN S A " V (VI. 1.6) 

The initial value of time is then defined before the loop which 
begins at statement 1. After augmenting the time, which is measured in 
hours by the variable TIME, this value is compared to 0.2 hr. That is, 
twelve minutes is considered time enough to accelerate to any practical 
value of Mach number. The lift coefficient is 

C = W/[.7 S„ P EM 2 ] (VI. 1.7) 

.L Wo 

where other input parameters are 

S ■ reference wing area, sq ft 
w 

W «* current weight, lbs 

This value of lift coefficient and the Mach number are input to 
subroutine DRAG, and the appropriate atmospheric properties are input 
to subroutine ENGINE. Together, these permit computation of thrust in 
lb, 
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T - P N EN p 


# 


(VI. 1.8) 


fu«l consumption in lb per hr, 


w 

FUEL 


w 

FUEL 



true airspeed in knots. 


(VI. 1.9) 


TAS - EM S /1.69 , 
A 


(VI. 1.10) 


and equivalent airspeed, 

EAS - TAS/SQRT (RHO/. 0023769) (VI. 1.11) 


where 

EN p ■ number of engines 

F ■ thrust per engine, lb 

N 

EM ■ Mach number 

RHO * a':, density at altitude, slugs per cu ft 
The acceleration is then expressed by Newton's Law as 

DVDT - (32 . 2/W) * (T - C *S *.7 P EM 2 ) (VI. 1.12) 

D W O 

where the drag coefficient C D is known from the earlier call to sub- 
routine DRAG. It is expected that the thrust exceeds the drag, so that 
this forwax'd acceleration is positive. If not, the subroutine merely 
writes a number of output quantities before returning. Otherwise, the 
velocity is augmented, 

• V - V + DVDT * DT (VI. 1.13) 


and the Mach number is 


EM - V/S R (VI. 1.14) 

A 

The aircraft weight and range are found, in lb and nm, as 

w - W - W puEL * DT/3600. (VI. 1.15) 

R - R + V * DT/6076. (VI. 1.16) 
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and tha fual consumed 


W p - W p ♦ Wpy^ DT/3600. (VI. 1.17) 

If tl current Mach number le leaa than the command value, tha 
time ia augmented at statement 1, and the loop la repeated. Otherwise, 
the final values of Mach number and true and equivalent airspeed are 
calculated corresponding to the input value, X.„, at statement 2, 
and a final line of output is printed. 
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V 1 . 1 . 3 Subroutine CUHB 

Thla Is a subroutine of 250 cards which la ^»i.lad In order to 

combine aerodynamic characteristics of the aircraft with Its thrust 

characteristics, to simulate the planar climb maneuver to the altitude 

^ENDCL’ The external subroutines CLIFT, DRAC, ENGINE and TP ALT are 

called, and three internal functions are defined in this calculation. 

Additional input quantities are constraints on pitch angle, altitude 

step sire, maximum airspeed and atmospheric properties. 

The subroutine first initializes a set of parameters, ending with 

an estimated value of maximum lift coefficient C The input parameter 

LL* 

ICLM specifies the type of climb maneuver desired, 

( 1, Climb at maximum rate of climb 

2, Climb at maximum allowable operating speeds (VI. 1.18) 
3, Climb at specified equivalent airspeed 


Following a writs statement, the loop which begins at Statement 1 
la the time-varying outer loop governing the climb maneuver. It begins 
with the call to TP ALT , which returns static pressure P o , gravity C 
and temperature T q corresponding to the current (initial) altitude, 

H (ft). The static pressure is expressed In lb per sq ft and the speed 
of sound is given in kts, 

S A - 49.1 (T ) 1/2 /1.689 (VI. 1.19) 

A O 

and the local air density in slugs per cu ft follows according to the 
equation of state, 

RHO - P /(S3..32 T G) (VI. 1.20) 

o o 

the minimum climb speed is expressed in kts and based on the estimated 
maximum lift coefficient, CLL, 


V r - 1.1 SQRT [W/ (1.426 S u RHO C)] 
L W LL 


(VI. 1.21) 
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where W is the current weight and S the wing area, and the factor 

W 

1.1 provides a margin above the stall speed. If this velocity is less 

than 50 kts, it is raised to 50. 

Following the setting to zero of sever;. 1 integer variables, the 

maximum allowable velocity at ' current altitude is calculated as 

limited by the allowable equivalent airspeed, V , the maximum operating 

Mach number, EM. . the speed limit below 10,000 feel altitude, or the 
MO 

cruise Mach number, 

The next lines of the subroutine related these constraints on the 
maximum velocity: 




V MAX " V/S 0 *T<RHO/.0023769) 

(VI. 1.22) 

(a) 

If 

3 CLM * 3 “ d V CLMB * V>’ 




V MM - V CUfB /SCRT(RHO/.0023769) 

(VI. 1.23) 

(b) 

If 

H < 10000. ft and V > 250 kts, 

MO 




V MAX " 250 -/SQRT(RHO/. 0023769) 

(VI. 1.24) 

(c) 

If 

V MAX ^ s a em M0 ^ s P eed corresponding to maximum 



Mach number) , 



where 


V - s EM 

MAX A MO 

The velocity V 2 is next defined as 


(VI. 1.25) 

! mln(V MAX* V CRU’ V LIMIT ) ' 

I CLM " * 

(VI. 1.26) 

( min ( V MAX* V CRU ) ’ 

X CLM " 2 ° r 3 


V CRU " SA EM CRU 


(VI. 1.27) 

V LIMIT * V limx /SQF.T(RHO/.0023769) 


(VI. 1.20) 
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and 

IEMcRU ■ cruise Mach number, from COMMON/UNIVX/ 

V, Tu _ • sea level limit speed, from COMMON /TOCLM/ 

LIMA 

All of these velocities correspond to steady-state horizontal flight, 
and are measured in kts. 

The loop begins with the dynamic pressure and lift coefficient, 

Q - 1.426 RHO V* (VI. 1.29) 

C L - W COS(RGAM 1 )/(Q* S w ) (VI. 1.30) 



Rate of Climb - V 2 sin(GAM) ■ V 2 (T - D)/W 


Figure VI. 1.2 Force Equilibrium in Steady Climb 


where the numerical factor accounte for the units of Vj (knots) and 
of Q (lb per eq ft). 
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The eubroutlne CLIFT return* the angle of attack, ALPHA, required 
for thia lift coefficient, and a rate of climb multiplier le defined 

as 

( 1. (DELTA U - 0) (VI. 1.31) 

PHI - { " 

l 1 + V (V -V )/ (11. 278 DELTAy) 


The parameter DELTA^ le initially at zero, but It later measures the 
altitude increment remaining to be climbed. The eubroutlne DRAG 
returns a value for Cp, and ENGINE returns a value for thrust per 
engine, T^, which is multiplied by the number of engines, EN p . The 
rate of climb le expressed In ft per min as (see Fig. vl.1.2). 


RC p0 W " 101,34 V 2 (T a' C D S W Q)/W (VI. 1.32) 

If the velocity V^ is too large, the rate of climb is negative so 
that the end of this iterative loop Is indicated by. 


V 2 - V 2 - 10. (VI. 1.33) 

RC 1 * RC P0W (VI. 1.34) 

If V 2 is now less than V^, as given by Eq. VI. 1.21, insufficient 
power is available for climbing, and the subroutine returns. Otherwise, 
another iteration of the rote of climb loop is performed at this lower 
speed, V 2< However, when RCp QM * s positive, control passes to 
statement 6, for computation of the flight path angle, 


R CAM2 " ARSINlRC poM /(101.34 V 2 ) 1 


(VI. 1.35) 


When this is within 0.1 deg of its value on the previous Iteration, 
convergence has occurred and control passes to statement 7. Similarly, 
the next seven lines are concerned with finding the flight path angle, 
R CAM4* at thc P* tc h attitude angle is a maximum at ^HEj^, 

and the corresponding rate of climb is HC 
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Assuming thst I CLM - 1, such thst the alrcrsft climbs st asslmum 
rate, the equations between statements 9 and 18 are devoted to an 
analysis of the rate of climb with the pitch angle constraint, ^tETA* 
relative to the rate of climb with the power constraint, RCp^y* 

This portion of the program involves the use of rather complex iterative 
procedures, including the three functions; 


part(x.x 1 .x 2 ,x 3 ,y) 


(X-X 1 )(X-X 2 ) 

"(X 3 -X 1 )(X 3 -X 2 ) 


(VI. 1.36) 


V l (U 2' U 3 HV 2 (U 3* U l )+V 3 (U r U 2L_ 

B02A(V ]L ,V 2 ,V 3 ,U 1 ,U 2 ,U 3 ) • 2 [V 1 (U 3 -U 2 )-*-V 2 (U l -U 3 )+V 3 (U 2 -U 1 ) 1 

aden(w 1 ,w 2 ,w 3 ,z 1 ,z 2 ,z 3 ) - z 1 (w 3 -w 2 )+z 2 (w 1 -w 3 )+z 3 (w 2 -w 1 ) 


(VI. 1.37) 

(VI. 1*38) 


When both of these constraints have been satisfied, control passes to 
statement 20, where the angle of attack is computed as 


ALPHA - THE MAJ( + EYEy - RGAM^ * 'RDTU 


(VI. 1.39) 


where 

EYE^ ■ input wing incidence angle, deg. 

The lift coefficient is then found by a call to CLIFT, and the drag 
coefficient follows from a call to DRAG. The total thrust required in 
the equilibrium climb condition is found as 


T a - W*PHI *RC TETA /<101.34 *V 2 ) +C d S M Q 


(VI. 1.40) 


and after division by the number of engines ENp, subroutine ENGINE 
returns the fuel flow rate per engine and other quantities related to 
the propulsion system. Following computation of the equivalent airspeed 


EAS, the altitude is compared with 
the integers 


H 


ENDCL* 


and until it is greater. 


"quo • h/del h * - 05 


(VI. 1.41) 


13 



and 


YALE 22 " N QUO + 1 


(VI. 1.42) 


are found, so that the altitude 

H 9 - YALE 22 * DEL h (VI, 3 

will be the next value of H, whether or not H^ la leas than 

The apparent purpone of thla logic is to develop an Integral nunber of 

additional computational stepa for the remaining altitude. 

The closing computations involves the routine computations for 


the rate of climb, 

RC - RC pow /PHI, (VI.l 

and the flight path angle, 

rcam 2 - ARS1N(RC/(101.34*V 2 )] , <vi -1 

which in degrees is, 

CAMMA - RCAM 2 * RTOD , {VI • 1 

and the lift coefficient, 

c L - w cos(rgam 2 )/(q*s w ) <VI>1 


The angle of attack ALPHA is then returned in degrees by a call 
to CLIFT, and the pitch angle is 

THETAp - CAMMA - EYEy + ALPHA . W 1 * 1 

Finally, the velocity is V ■ V^. 


.43) 


.44) 


.45) 


.46) 


• 47) 


.48) 
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VI. 1.4 Subroutine DLAND 


This subroutine models the landing maneuver in four segments; a 

glide frott 50 ft altitude, a flare to touchdown, a delay before braking 

and a braked roll. Figure VI. 1.3 shows the flare and touchdown schematics. 

The program begins with the specification of landing weight, which 

is determined by the input integer I... . which is defined as 

WLD 


*WLD 


0, landing weight equals gross weight 

1, landing weight equals weight at end of previous 
mission segment 

2, landing weight equals input percent of gross weight 



Figure VI. 1.3 Flare and Touchdown Maneuver 
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The next five cards specify the landing weight (lb) and the wing 

loading WOS (lb per sq ft), and the initial altitude is 


H 


IN 


APP 


♦ ALT 


LND 


(VI. 1.49) 


where 

H^PP ■ approach altitude above runway, ft 
ALT^p ■ altitude above sea level of runway, ft 

The subroutine TP ALT again provides static temperature T q and pressure 

p at the initial altitude, which may then be modified for an increment 
o 

TDEL to the temperature, and by converting the pressure to lb per sq ft. 
The air density in slugs per cu ft is 


RHO = P / ( 53 . 32 G T ) 
o AMB 

where the ambient temperature, in degrees Rankine, is 


(VI. 1.50) 


and 


AMB 


T ♦ TDEL 
o LD 


G = acceleration of gravity, 


ft 


2 

per sec 


(VI. 1.51) 


The square root of the density ratio is useful in computing the equiva- 
lent airspeed at sea level. This factor is denoted by 

SR DRAT “ SQRT (RHO/. 002378) (VI. 1.52) 

Next, the flight path angle RGAM^, the angle of attack ALPHA, 
the approach altitude In 'wing spans HOB, and the stall Mach number 
EM are defined, prior to a call of subroutine CLIFT, which returns 
the appropriate maximum lift coefficient The associated stall 

speed in ft per sec is given by 


V gT - SQRT (2. WOS/ (RHO * C ) ) 


(VI. 1.53) 
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and this is converted to the equivalent airspeed 


V STEAS " SR t>RAT V ST 


(VI. 1.54) 


If the aircraft is single engine, and this velocity is less than 
the maximum stall speed of 61 kts or 103.03 fps, per the Federal Air 
Regulations Part 23, the maximum wing loading is calculated and printed 
as 


" CS MAX * 12 ' 6212 C UB 


(VI. 1.55) 


where the numerical factor is the dynamic pressure at sea level, In 
lb per sq ft, corresponding to this stall speed. 

The glide velocity in ft/sec is next defined as the product 


V ■ V V 
GL RAT ST 


(VI. 1.56) 


where 


V ■ input approach speed margin (default value » 1.3) 


and the height above the ground in wing spans at the beginning of the 
glide is expressed as 


H 


OB 


- W» 


(VT.1.57) 


where 

H^pp * approach altitude, ft 
B ■ wing span, ft 


The Mach number and Reynolds number per unit length are given by 


m a - V ct/ S ‘ )RT < I - 4 * 0 * 53 ' 32 * T AMB ) 


(VI. 1.58) 


where 


RE 


LI 




■ kinematic viscosity, as returned by TP ALT. 


(VI. 1.59) 
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Next, the eubroutlnc AERO le celled, which returns e number of coeffi- 
cients epproprlete to the dreg coefficient of the aircraft of Interest. 
The dynamic pressure is calculated by the usual equation, and If 
I ENCS T " subroutine ENGINE returns a value for Idling thrust, 

The flight path angle iteration loop begins at statement 7, where 
the lift coefficient is computed in terms of the flight path angle, 


<* L " WOS * COS(R ( j AM 1 )/Q ( jk 


(VI. 1.60) 


Here, R CAM1 ia the variable being calculated In this loop, and It 
is positive when the sink rate is positive. The wing leading Is WOS 
lb per sq ft and the dynamic pressure corresponds to the velocity, 

Vql* The drag coefficient appropriate to the approach configuration 
la returned by subroutine DRAG, and the positive rate of sink is 
expressed for the idle thrust as 


RS 

IDLE 


60.0 V 


GL 


(C S Q - T 
D W W GL 


IDLE 


EN p )/W L 


(VI. 1.61) 


where the numerical factor is to get the units of ft/min. 

This equation is followed by another for the rate of sink RS, in which 
the actual thrust is expressed as T . The equilibrium flight path angle 

r\ 

associated with RS is given by 


RGAM - ARSIN [RS/60. 0 * V ] 
2 GL 


(VI. 1.62) 


and if this is within .1 deg of the previous value, convergence has 

occurred. Otherwise, RGAM^ * RGAM^ and control returns to statement 

7, where the lift coefficient is recomputed according to Eq. Vl.1.60. 

After this loop has converged, several constraints must be checked, 

which may require the loop to be re-entered. The sink speed is checked 

to see if it is limited by the idle thrust value found in Eq. VI. 1.61 or 

exceeds the maximum allowable value RS „ . If the maximum allowable value 

MX 

is exceeded, the necessary thrust increase is computed to provide the 
desired sink rate, 


T A " °D S W C GL - \ ^m/ 60 ' 0 V GL ) 


(VI. 1.63) 
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before reconverging to the flight path angle RGAM^ • 

At statement 15, the flare maneuver is Initiated with computation 
of the initial pitch angle and the final flight path angle at touchdown, 

THETA ■ ARS IN IRS/60. 0 V )) (VI. 1.64) 

GL 

°"* TD • SI V a ' (VI.1.65) 


where 


SINK TD " vertical velocity at touchdown (3 ft per sec or 
input) 

The input variable XLF indicates the maximum allowable load factor 

MX 

in g's if less than 4., and the radius during the flare is computed as 


» 2 ■ v 0 l /(g(xu, mx - l -» 


(VI. 1.66) 


and the altitude at which the flare is initiated is approximated as 


h flar - v gl <TH m2 - g “to ,/ 12 - - l -" 


(VI. 1.67) 


Otherwise, the flare radius is 


R - V / (G (XLF - l.)J 

I* bL 


(VI. 1.68) 


where 


XLF - V^ L (THETA 2 -GAM^ D )/(2*C*H FLAR ) 4 1. 


(VI. 1.69) 


However, if the input variable X1J\,„ is greater or equal to 4 it 

MX 

is taken as the flare height in ft. 


Hflar 


XLF, 


(VI.1.J0) 


MX 


The computation continues with the air distance from approach 
altitude H App to intersection with the ground, if less than 

h app' 
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(VI. 1.71) 


DL CL • "appear (THETA) 

and from this point to touchdown (••# Fig. VI. 1.3). 

DL tr - R z THETA/2 . (1 - CAM^ /THETA) 2 (VI . 1 . 72) 


The latter expression is a small-angle approximation to 


DL tr • R z |sin(THETA)-sin(CAH TD ) + 


cos (THETA) -cos 


L 


tan (THETA) 


m 


If **plar ^ceeds H^p p , DLg L * 0 and the air track distance i», 
to first order, 

DL tr - K,(SQRT(2 H app /R 2+ CAM^ d ) -CRH^I <«•!. 

The touchdown velocity is expressed as the average of the stall 
speed and the approach speed 


^TD " (V RAT* 1#)V ST /2, 


(VI. 1.74) 


where V ) 1 is the ratio of approach speed to stall speed. The 

KA I 

delay distance, over which the velocity changes, is 


D - V T 

DELAY TD DELAY 


(VI. 1.75) 


where T DELA y is input or has the default value 1 sec. 

Other parameters which are then computed are the height in wing spans 
above the ground and the Mach number at touchdown. These are given 
respectively by 


HOB - H /B 
TG 


(VI. 1.76) 


and 


EM l “ V td /SQRT(1.4 * G * 53.32 * T^) 


(VI. 1.77) 
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whore 


H - altitude at touchdown 

TG 

These are required Inputs to the subroutines CLIFT and DRAG, 
which are next called. They return the gear-down lift and drag coeffi- 
cients, C LRL and The net drag is the sun of aerodynanlc drag 

and friction drag, which Is proportional to the difference of weight 
and lift. A one-dimensional equation of notion takes the fora, 

5t " v 5* * w * T ' D ' ,1<M " w * 


where the symbols have standard definitions and where |i Is the coef- 
ficient of friction. In the present case, D and L are both propor- 
2 

ti~nal to V , and the equation becomes separable, as follows: 


VdV 


F RAT* A RAT v 


2 " ~ 2 ’ A RAT dX 


where 


A RAT " 

dl rl /ic lmx v st 1 

(VI. 1.78) 

F ■ 

MU_ - TOW. 

(VI. 1.79) 

RAT 

B L 




(VI. 1.80) 

■ 

es 

Q 

C DRL “ ^B C LRL 



and 

MU ■ friction coefficient 
B 

TOW » thrust to weight ratio at landing 
L 

The equation can then be integrated, to give the ground roll 
distance X as 


r Frat 1 

„l 0 - -n.o2e7[w os /(o MT V> 108 [viJ 


(VI. 1.81) 
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where the numerical factor la tha invar a* of the product „ «a level 
air density and gravity. The deceleration during the ground roll 
maneuver in g'a is the quantity 


OG 


V «' (J 


* G * 


D LG> 


Finally, the total landing diatance ia the sum. 


(VI. 1.02) 


DL - DL_ + DL ♦ D_ „ ♦ DL_. 
T G TR DELAY GL 


(VI. 1.83) 


The remaining atepa in the subroutine relate to satisfying 
a runway length constraint by iterating on the airplane's wing 
loading, w . This logic is dependent on the input value of runway 

viD 

length requirement, X T ___ . and it requires two iterations on the wing 

laDGRQ 

loading. An interpolation or extrapolation of the associated runway 

lengths then develops a required maximum wing loading. 

The procedure is described as follows. Initially I » 1 and 

LER 

I » 1, when DLAND is first called by MAIN. If DL is found to be 
WGS T 

less than the required value of runway length, the only remaining compu- 
tations deal with conversion of velocities to knots and with augmenting 
the runway length to meet FAR Part 25/121 regulations; i.e., the required 
runway length in ft is 


D 


FAR 


DL t /.6 . 


(VI. 1.64) 


On the other hand, if exceeds the limiting value, x LDGR q, 

the initial values of wing loading and runway length are relabelled W cgp 
and X..., and the wing loading is reduced by 30%; 



.7 



(VI. 1.85) 


After augmenting I WGS , and setting th-s subroutine argument to 


I LER 


2 


(VI. 1.86) 


control is returned to the MAIN calling program. On the following itera- 
tion, I^g - 2 and the shorter tunway length DL^ is found for the 
reduced wing loading, V cs . The sensitivity is found as 
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SLOPE • (W - W }/(X - DL ) '**' 

1 GSP Gs"' PAX T # 

and the counter I tl __ is sunmentad to 3, while l.__ ■ 2. This slope 

WvjS Lla 

is positive unless an erre* has occurred, and the required wing loading 
is found at statement 75 as 


W 


GS 


- W Gsp ♦ SLOPE * 


iX LDGRQ " X PAS ) 


(vi.; 


Finally, if the runway length requirement is met, no further 

iterations are needed, and at statement 100, the counters 1 and 

LER 

*WGS are 8et to ^ ' and the concluding computations ending with Eq. 

VI. 1.83 ore performed. 


. 87 ) 


. 88 ) 
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Subroutine TAKOFF 


VX.l.i* 


The simulation of the take-off maneuver requires consideration of 
several smaller segments or '"events." Each event is terminated when 
a particular performance variable reaches a given value. These events 
are summarized as follows, with key variables in parentheses: 

a11 engine j at take-of.t power: 

(1) Begin ground roll (velocity) 

(2) Rotation (velocity) 

(3) Lift-off (lift * weight) 

(4) Begin : ar retraction (altitude) 

(5) Altitude - 35 ft 

(6) Altitude - 50 ft 

(7) Begin flap retraction (altitude) 

(8) End of take-off (altitude) 

One engine out: accelerate-stop distance 

(1) Begin ground roll (velocity) 

(2) Engine failure (velocity) 

(3) Remove power (time) 

(4) Apply brakes (time) 

(5) Stop (velocity) 

The take-off computation is represented by simultaneous integration 
of the differential equations for velocity, flight path angle, distance 
and altitude, using a time interval of 0.2 sec. 

The subroutine DERIV, to be discussed in the next subsection, per- 
forms the computation of the time derivatives, and at each time step, 
updated values for weight, thrust, air density, etc., are found before 
the integration subroutine INTS is called. ^ 

Subroutine TAKOFF begins with the definition of a linear interpola- 
tion function of the form. 


YYY(X) - Y x 




V x i 


<w 


(VI. 1.89) 


which defines a value YYY between Yj^ and Y 2 » when X falls between 

X and X . This is followed by initialization of some 30 numerical 
12 

factors, ending with the vector values T(l)-T(ll), used by the integra- 
tion routine 
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The subroutine TPALT then returns staospherlc properties st the 
Input altitude of the airport and runway, and the local speed 

of sound and air density are given as 


S A - 49.1 SQRT(T amb ) 

RHO - P o /[53.32 G] 


where 

T,,_ ■ ambient airport temperature, deg Rankine 

AHo 

2 

G - local gravity, ft per aec 

P o - local static pressure, lb per sq ft 

The sea level stall speed is given in kts by 


V STLKT 


SQRT 


2W 

•° 023769 S W C LMAX 



(VI. 1.90) 
(VI. 1.91) 


(VI. 1.92) 


where the symbols have their traditional definitions. 

c - c A 2 

LTO LMAX RAT 


(VI. 1.93) 


is determined where v is the input takeoff speed margin or has a 

KAi 

default value 1.10. 


Speed increments are then used to define the engine failure speed, vi, and 
the rotation speed, VR, 


V + 

STLKT . 


DV. 


(VI. 1.94) 


V R - v i + DV r (VI. 1,95) 

where 

DV^ * increment of "decision" speed above stall, kts 

DV ■ increment of rotation speed above V , kts 
R 1 

The loop between statements 180 and 186 finds the speed VEL for the 
best rate of climb, beginning with the stall velocity. This loop uses 
lift coefficient C and the Mach number EM, as inputs to subroutine DRAG, 
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which returns the drag coefficient, Cp. Subroutine ENGINE returns 
the thrust available, T , and together these permit computation of 
the power required, 


pw v 

and of the power available, 


5 RHO S,, VEL c * VEL 
W D 


pwr avl - t a en p vel 


(VI. 1.96) 


(VI. 1.97) 


The velocity is augmented by DEL^ - 10. feet per second, and the loop is 
repeated, until the difference PWR AVL - PWR r ^ r is less than it was on the 
previous iteration. The velocity increment is then reduced to 1 foot per 
second, and the velocity is eventually found at which the available power 
exceeds the required power by the maximum amount. This permits computa- 
tion of the best rate of climb velocity in kts, 


V END " * 592087 VEL^ * SQRT (RHO/. 002 3769) 


(VI. 1.98) 


which must satisfy the constraints, 


V END * 250. kts 


(VI. 1.99) 


and 


V < V 
END - CLMB 


(I cu. • » 


(VI. 1.100) 


The next important portion of the subroutine deals with the loop 
between statements 1 and 300, which treats the ground roll, rotation, 
lift-off and climb phasep of the take-off maneuver. The sequence begins 
with the ground roll equation of motion, for which control passes to 
statement 230. The Mach number and altitude are input to ENGINE, which 
returns the thrust and fuel flow. The horizontal acceleration is the . 
expressed as 

T(8) - (32.16/W)(T a cos(ALPHA r ) - D 


- MU * (W - L - T a sin(ALPHA R ) ] 


(VI. 1.101) 
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where all of the terras are tine varying during this phase, except the 
angle of attack and the friction coefficient. Control passes back to 
statement 225, after setting 1^ ■ 1, and the ground roll Integration 
continues after a call to INTM. Until the velocity exceeds V R , the 
Integration proceeds with T(5) ■ flight path angle • 0. 

When the rotation speed is exceeded, Interpolation develops values 
for the time and value of the lift-off speed, TIMR and VW, using the 
function YYY defined in Eq. VI. 1.89. This requires setting I R • 1, 
and the lift and drag coefficients C^ R and are also found at this 

time by linear Interpolation. The time after the rotation, when a more 
or less steady positive flight path angle has been acquired, Is 

TIM^ - TI^ + DELT yR (VI.1.1P2) 

where DELT is input, or defaulted to 3.5 sec. During this brief 
V R 

interval, of course, the exact dynamic response of the aircraft would 
require consideration of the aerodynamic derivatives, the ground effect 
on these derivatives and the aircraft inertia in pitch. These details 
are not considered in this subroutine, which is intended to model the 
aircraft as a jxjint mass during the take-off maneuver. 

Following the computation of angle of attack, it is tested whether 
the vertical forces (lift and upward component of thrust) exceed the 
weight, with T(5) - 0. If so, the friction is set to MU ■ 0, and 
the take-off speed is calculated in kts as 


V T0KTS ’ T(4) * * 592087 (VI. 1.1 

The parameter values at lift-off are Initial conditions for the fourth 
event (gear retraction) and this is modeled by a reduction in drag coef- 
ficient. This begins when the altitude equals H^ R “ 20 ft, and the 
corresponding time is found by Interpolation as Tg. The gear retracting 
operation requires an Interval of DT (;il ■ 7 sec. and the gear drag incre- 
ment falls linearly to zero over this time Interval. The flag is 

set to 1 at the beginning of retraction, and the lift varies linearly 
with time from C LR to C LT(J , at beginning and end of gear retraction. 
This permits computation of angle of attack during gear retraction. 
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The next two events correspond to take-off altitude heights of 
35 ft and 50 ft, and these occur at runway distances found by interpola- 
tion as and ft, when I pAS ■ 1 or 2, respectively. When this 
flag has been set at 2, the altitude exceeds 50 ft and control passes 
over the 23 statements dealing with these two events. 

The flaps are deflected at the angle DFLP , and are to be re- 

TO 

tracted at the input altitude H FLp * 400 ft. Retraction requires a 

time interval which is calculated as DT , . which is defined as equiva- 

lent to a retraction rate of 3-1/3 deg per sec. Over this time interval 

lift and drag characteristics vary linealy between flaps down and flaps 

up conditions. The take-off maneuver terminates at the specified 

altitude H ... If H.,,„ is less than H„, „ it terminates when the flaps 
MAX MAX FLP 

are retracted. 

Engine failure is simulated by repeating the ground run portion of 
the subroutine, when the input flag is N ■ 0. 'This failure is 
assumed to occur at the ground speed, V^, where all parameters are 
calculated by interpolation. The flag I M is set to 2, indicating engine 
failure has been initiated, unless the total number of engines is 
EN p = 1. The continued takeoff computations after engine failure is 
handled by the same equations as the all engine takeoff except the 
thrust is reduced by the loss of the engine. The continued takeoff 
calculations are terminated at an altitude 50 feet above the airport. 

The accelerat-stop computations after engine failure begin at 
statement 400. Following initialization, the loop beginning at state- 
ment 920 integrates the equations of motion until the brakes are 
applied or until the velocity is less than 1 ft/sec. This is done by 
implementing the following logic: 

(1) Engine thrust is reduced to zero 1.8 sec after failure is 
initiated. 

(2) Brakes are applied 3.5 sec after failure. 

A closed form solution of the one-dimensional motion is then used 
with the braking value of friction, MU , and this computation has been 

B 

described in subroutine DLAND. The result is the incremental runway 
distance DL_, which is the distance traveled after brakes are applied. 
The total runway distance is then given by 
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T (6) ♦ DL 
G 


(VI. 1.104) 


D 


STOP 


# 


where T(6) is the runway distance preceding brake application. 
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VI. 1.6 Subroutine DEFIV 


Thia 90-card aubroutina la called by TAKOFT, aa an argument of sub- 
routine I NTS, and It computes the time derivatives of the aircraft 
position and velocity. The two kinematic equations are Integrated to give 
horizontal and vertical distance traveled, and the two acceleration equa- 
tions in general yield the velocity and the flight path angle, while 
Ignoring the rigid body pitch dynamics of the aircraft. The acceleration 
and velocities are limited by input values of the following quantities: 

(1) Maximum load factor 

(2) Maximum velocity 

(3) Maximum lift coefficient (angle of attack) 

(4) Maximum pitch angle 

(5) Minimum (zero) load factor 

(6) No de.-oceleration along flight path 

In simplest form, the forces acting on the center of mass of the 
aircraft are lift, drag, thrust and weight, and if the thrust is oriented 

by the angle o - 1^, the acceleration components acting on the aircraft 
mass are 


L 



Figure VI. 1.4 Applied Forces During Take-off 


30 


(VI. 1.105) 


0 - (g/W) (T cos (a - 1^) - D - W siny - 

V y - (g/W) [T sin (a - iy) L - W cos y ) 

where M ■ 0 when the aircraft is airborne, and where the normal force 
is otherwise: 

F - W - L - T sin (a - i ) . (VI. 1.106) 

N W 

The subroutine accounts for weight loss due to fuel consumption , and the 
thrust varies according to the speed and altitude, as found by sub- 
routine ENGINE. The angle of attack is equal to the wing incidence 
angle until the aircraft rotation is initiated during takeoff. 

The subroutine begins with the modification of the angle of attack, 
ALPHA, during the airborne portion of the take-off maneuver. As measured 
by the load factor the angle of attack lu degrees Is Increased 

according to: 

ALPHA + .1 \f > * 85 

ALPHA + .3 .8 < X^ F < .85 

ALPHA ♦ .5 ~ ' 8 

Whether or not the aircraft Is airborne, subroutines CLIFT and DRAG 
are called In mode 2, returning values of and Cp. The next 12 lines 
of code deal with the drag and lift Increments due to landing gear and 
flaps. These events occur at the times (begin landing gear retrac- 

tion) and Tpjj, (begin flap retraction), and quantities returned by 
DRAG are 


ALPHA - 


Cpcu " drag coefficient, takeoff flaps down and gear up 
C DFU " ^ rag coefficient, flaps up and gear down 


The corrected coefficients of drag and lift are then 


C D ’ C D ^DG G FAC /DT GR " DC DF F FAC/DT, 


FLP 


C ■ C - DC F /DT 
L L LF FAC FLP 


(VI. 1.107) 
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where DC^ and DC^ are increments in drag due to gear and flaps and 
DC is the lift increment due to flaps. The remaining factors repre- 

LF 

sent the linear change of these increments over the time interval 

DT „ or DT . 

GR FI.P 

Maximum and limit values of lift coefficient are then found, to 

account for takeoff speed margins. If C exceeds the limit value 

L 

C . the subroutine CLIFT returns a lower angle of attack, and a 
LLMT 

second iteration from statement 1010 develops a lift coefficient less 
than 

The fuselage actltude angle may also limit the angle of attack. If 
THETAp exceeds THE^^. In this case, the angle of attack is defined 
in degrees as 


ALPHA - THE^ + EYEy - T(5) * 57.29578 


(VI. 1. 


where T(5) is the flight path angle. Control again returns to the 
beginning of this loop, at statement 1000. Otherwise, the normal load 
factor Increment due to lift and thrust is found as 


*LF " [C L * QSW + t A sin(ALPHA R )J/W cos(T(5)) 


(VI. 1. 


and if this exceeds when the aircraft is airborne, the angle 

of attack is reduced by .1 deg for the next Iteration. This reduction 
is only .05 deg If the flaps have already begun to be retracted. 

The equivalent airspeed is then calculated in kts as 


EAS - T(4) * KTFPS * SQRT(RH0/RH0 $L ) 


(VI. 1. 


where 


KTFPS - .592087 and RHO 


SL 


.0023769. 


This is required to be less than the velocity for maximum rate of 
climb, V END * Otherwise, T(4) is reduced to this equivalent value and the 
velocity acceleration is set to T(8) ■ 0. 

The tangential acceleration in the unconstrained case is then given 
in statement 4200 as 


108) 


109) 


110 ) 
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(VI. 1.111) 


T<8) - (32.15/W) (T co«(MJ>HA o ; - D - MU * P - M*SIN(T(5) ) ) 

A R 

vbari drag and noreal forca during takaoff ara 

D - C D * QSW 

p • W - * QSW - r. * SIN(ALPHA r ) (vx. 1.112) 

If » 2 or T(8) le positive, the normal acceleration T(9), 

and horizontal and vertical velocity components T(10), T(ll), are 
computed by straightforward aquations before returning. If T(8) Is 
negative, the drag exceeds the thrust, so the angle of attack Is re- 
duced and the counter 18 augmented. If 20 Iterations do not 

develop a positive acceleration T(8), a lower flap deflection is sug- 
gested by the print statement when I ALQp " 20. 
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VI . 1 . 7 Subroutine TAXI 

This is a very brief subroutine with the purpose of computing the 
fuel burned while taxiing at ground idle. This is used to account for 
taxi time at the start and end of mission. This fuel weight is sub* 
tracted from the initial gross weight to provide the aircraft weight 
at take-off. 

At each time step, the subroutine is called, and it calls TPALT 

for values of pressure and temperature, P.„_, T._,_, at the runway 

AMB AMB 

altitude, H. These are input to the subroutine ENGINE, which returns 
fuel flow per engine, FF, in lb per hr. This is multiplied by the 
number of engines, EN p , 

FF - EN p * FF (VI. 1 

Next, the weight of the fuel burned during the time interval DELT^ 
is 

WF T - FF * DELT t (VI. 1 

and this reduces the aircraft weight, 

W - W - WF_ (VI. 1 

T 

The total fuel burned and the total taxi time are then augmented as 

WF - WF + WF T * (VI. 1 

and 

ST - ST + DELT_ (VI. 1 

T 

If Kyujjg is nonzero, both input and outpue values of these time- 
varying parameters are printed before control is returned to subroutine 
PERFRM. 


. 113 ) 

. 114 ) 

. 115 ) 

. 116 ) 

. 117 ) 
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VI. 1.8 Subroutine TORN 


This la another very brief subroutine, which calculates certain 

parameters descriptive of the steady turn naneuver. The subroutine 

begins by determining the static pressure P and temperatura T 

o o 

at the input altitude by calling subroutine TP ALT. The speed 

of sound is then given in ft per sec by 


A SON * 49 - A SQRT( V 


(VX. 1.118) 


and the dynamic pressure in lb per aq ft is 

- .7 * P e * EdJ (VI.1.11 

where EM^ Is the airplane's Mach number in the turn. 

The load factor in the turn la the ratio of lift to weight, as shown 

in Figure VI. 1.5. With the load factor XLF TRN input as greater than 
one, this permits the required lift coefficient to be found as 



Figure VI. 1.5 Aircraft in Coordinated Turn 
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(VI. 1.120) 


Saqt * ^TRN ™'%m * 


where WOS le the wing loading In lb per sq ft, at calculated a few 
llnee earlier. If thla la leea than the Unit value 
angle le found. Otherwiae, the lift coefficient is C 
limit load factor la 


C LTLMT* th * bank 


'LTLMT 


and the 


^TRN " C LRQT ^TURN^ W0S 
In either caee, the bank angle la given by 


(VI. 1.121) 


PHI - ARC0S(1/XLF TRN ) (VI. 1.122) 

as suggested by Figure Vi. 1.5. The Mach number and lift coefficient together 
Imply a drag coefficient c DTU rii* aa provided by subroutine DRAG, and 
the minimum thrust required for maintaining velocity during the turn is 
equal to the drag, 


TR TRN " C DTURN Q TURN S W (VI. 1.123) 

T he steady turn radius In ft la then found by equating horizontal 
forces to the centripetal acceleration, or 

tiu, iiad ■ a son /[0 ‘ xlf ibk SIN(PHI)] m.i.iu) 

Finally, If I DC is not equal to 99., subroutine ENGINE determines 
the thrust available at this flight condition, and if It Is less than 
that required, an appropriate message le printed before returning to 
subroutine PERFRM. 

3 f I * 99 subroutine TURN is being used by subroutine ENGSZ to 

DC. 

determine the thrust required in sizing the engines to the specified 
turn requirements. 
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VI. 1.9 Subroutine 1 XRANGF 


Title ■ubroutlne computes the cruise range using the Brequet range 
equation, for both propeller and turlcjet aircraft. The equilibrium* 
in cruise, of lift and weight and of thrust and drag together lead to 
cruise angle of attack and cruise velocity, for a particular power set- 
ting. The equilibrium equations are nonlinear and must be solved itera- 
tively. and subject to several constraints. 

The subroutine begins by determining static pressure P Q and the speed 
of sound A at the input cruise altitude, H, through a call to TP ALT. 
Initial values for slrcraft weight and Mach number are also specified, 
and the square root of the density ratio Is defined as 


SQRDR 


[P /P 
o 

T /T 
. o 


ll/2 


REF 

REFj 


(VI. 1.125) 


where Pgpp " 2116.22 lb per sq ft and Tgg p ■ 518.67 deg R. This scale 
factor will be needed at several points la'ter in the program. 

The subroutine ASPEED is called, the engine operating condition 
for normal cruise power having been specified with the indicator K EN q* 

The Mach number returned as EM^ is the "normal power" value, and it 
corresponds to the velocity in kts, 


VEL r - 29.07 EM^jp SQRT(T Q ) 


(VI. 1.126) 


This also serves as an initia). value for the loop between statements 20 and 
30, which leads to the speed for the best "specific range," in nautical 

miles per lb of fuel; 

SPEC R - VEL K /(FF EN p ) (VI. 1.127) 

• 

where FF is the fuel flow returned by subroutine ENGINE, and EN p is 
the number of engines. This iterative loop varies the velocity, first 
in steps of 10 kts and, once the peak is surpassed, in 2 kt step* . The 
resulting "best range" Mach number is denoted . 


t 
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The next phase of the subroutine deals with cruise performance at 
a maximum of three fuel capacities and at three Mach numbers. The fuel 
weights are input through COMMON /WGTRCE/ as 

W FAMAX " oaximuo available fuel weight (minimum payload) 

W FAMIN " minimu ® available fuel weight (maximum payload) 

W FADES « design available fuel weight 

The Mach numbers have been specified as 


£M x 

em brg 


cruise Mach number. Input as cr j^qj , or defaulted to EM CWJ 

normal power Mach number 

maximum specific range Mach number 


The parameter SPEED used by subroutine ACCEL is equated to one of 
these Mach numbers, depending on the value of the input flag and a 

number of other definitions are made at this point. 

The subroutine ACCEL returns the time ST, the range P, and the 
fuel weight DEL^^ required to accelerate from the "climb" Mach number 
EM C lb t0 the "cruise" Mach number SPEED. The integers I DC and Ig EG ^ 
are defined as 1 an.^ 9 for later use in the subroutine *'1GINE, and 
weight of available fuel is defined as w p^vLB’ w *’ ich depends on the 
mission weight and on W FAM ^ X * W pAMIN and W FADES* The ' ,mls8ion weight," 
^MISN’ lf input is a take-off weight which is less than the gross weight 
but greater than the empty weight of the aircraft. 

After setting to zero the parameters OUT(I,J) and RES (I) , the 

F 

index K is initialized at zero, before the loop 120 is entered. This 
loop is passed through for each of the three Mach numbers, EM (K=l) , 

A 

EM^p (K=2) , and EM br( , (K=3) . Initially, the Mach numbers are checked 
to determine if they are within /he operating envelope of the aircraft. 

The velocity V is expressed in kts, and if the operating velocity 

MAA 

V exceeds V , then V v and EM are re-defined to correspond to V.... , 

X MAX X MAX 

which must satisfy the maximum operating airspeed and Mach number 
constraints. 
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For all but fixed-pitch propeller eye terns (i.e., for (CODE 4 8), sub- 
routines CTAER and ENGINE are called to determine the drag DRG and the 
required thrust of the engines, FN_ n , at the flight condition of inter- 

KVJ 

eat. These parameters are defined at statement 145 at the altitude, H, 
and aircraft weight, W. 

A longer, iterative process is required for fixed pitch propeller 
systems, or KODE ■ 8. In this case the aircraft speed is determined by 
the engine *>ower setting, to find a Mach number SMN at which thrust and 
drag are in balance, i.e., at which 


ERROR - (FN - DRG) /DRG < .002 , 


(VI. 1.128) 


where FN _ n is the thrust as output by ENGINE and DRG is the drag as 
output by CTAER. 

On convergence of this 2oop, if shaft horsepower SHP is less than 
available horsepower (HP^^g) , control passes to statement 146. Other- 
wise the engine RPM Is reduced by 52 and the loop is reinitiated. This 
can be repeated until convergence occurs. The equation.! between state- 
ments 146 and 145 are concerned only with the time spent in accelerating 
to the cruise Mach number, and with the fuel consumption during this 
interval. When KODE - 8 and 1^ ■ 0, as on the initial arrival at 146, 
subroutine ACCEL returns the aircraft weight and range covered during the 
constant altitude acceleration to the Mach number SPEED. This is necessary 
regardless of the Iterations needed in the Initial loop between state- 
ments 142 and 146. 

The Breguet range constant at the start of cruise is then found in 
units of nautical miles, as 


BREG 1 - V x W x 3600. /(FF X * 6076.1) 


(VI. 1.129) 


where aircraft velocity in ft per sec is the product of Mach number and 
the speed of sound. 
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(VI. 1.130) 


V x “ S KN A S0N • 

aircraft weight la lb and fuel flow la, 

FF X - FF * EN p (VI. 1.131) 

where FF la fuel flow per engine in lb per hr, and the aubacrlpt 1 
Indicates parametric valuea at the beginning of the cmlae phase. 

The fuel weight available for cruise la the difference 


where 


w ■ w •> Vf « U 
FCRU FAVLB FI FRES 


(VI. 1.132) 


W pl * fuel weight used in taxi, takeoff, and climb, lb 
W - reserve fuel weight, lb 

r Khb 


Aircraft and fuel weights at the end of cruise are Wj and W p 2 » 
respectively, and these are found by the same Iterative procedure as 
described above for the beginning of cruise. . The subscript 2 indicates end 

of cruise and is used to specify fuel flow FF 2> Mach number EM^ , velocity. 


V 

KTS2 


.5921 EM 2 


a son 


(VI. 1.133) 


and equivalent airspeed 


V KES2 ’ V KTS2 SQRDR ' (VI.1.134) 

where 8QRDR is defined in Eq. VI. 1.125. 

This permits the computation of the end of cruise Braguet factor 
BRKC. by an equation similar to that of Eq. VI. 1.129, following which the 
"average" factor is 


breg a - (breg jl ♦ brec 2 )/2 


(VI. 1.135) 
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Th# cruise rang* la than ealculatad aa 


R cww ■ breo a awcii./U-Vm/VI 

and other parameters at the end of cruise are given as range R 2 , 
angle of attack ALPH 2> attitude angle THET^* lift coefficient C L2 , 
lift to drag ratio ZLOD 2 , and time STj. 

The array 0UT(6,16) Is then used for storing computed quantities 
at start and end of cruise for the three Mach numbers* as Indicated by 
the counter K. A large number of WRITE statements follow, which are 
self-explanatory, and these are followed by a redefinition of available 
fuel; i.e., the first pass through the loop beginning at line 100 is 
with the available fuel weight given by 


W FAVLB " W FAMIN " DEL WGT 


(VI.l 


while the subsequent passes define this weight as W_. UAV 

FArnX 

respectively. 


and 


W FADES ’ 


The closing computations deal with the specification of the output 
parameters, at the design payload fuel weight, computed on the final 
pass through XRANGE, which are used for checking range or endurance 
requirements in subroutine RGBAL or operating cost in subroutine GACOST. 

The input parameter I CRUS specifies the speed condition selected 
for the design mission performance and the index 


KK = 


2(1 + 
CRUS 


1 ) 


1 , 


and the following 


R 


CLBNM 


R. 


CRNM 


CRHR 


CLHR 

em save 


parameters are defined: 

range after climbing, to start of cruise, nm 

cruise range, nm 

time spent in cruise, hr 

time spent in climbing to start of cruise, hr 
Mach number at start of cruise 


.136) 


.137) 
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"save 

W FCRU 

R 

ST 


■ altitude of cruloo, ft 

■ fuel consumed in cruise* lb 

■ total range* after cruise* tun 

■ total time from take-off* after cruise* hr 

■ total fuel used, after cruise, lb 
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VI . 1 . 10 Subroutine ASPEED 


This subroutine is used to calculate the maximum Mach number for 
equilibrium flight at given values of weight, altitude and engine power 
setting. 

The subroutine is called by subroutines GACOST , X RANGE and PNOYS , 
and it calls on subroutines CTAER, ENGINE, ITRMHW and TP ALT. The itera- 
tive method begins with an estimate of Mach number, to which corresponds 
a thrust value based on the engine data, at a specified engine power 
condition, KENGX. Subroutine CTAER then determines the cruise angle of 
attack at which lift and weight are equal, and returns a value for drag, 
which is compared with the thrust before modifying the Mach number. 

The call arguments of ASPEED are the output Mach number XM , alti- 

N 

tude (ft) H in , weight (lb) W IN > and engine power indicator After 

specification of certain parameters, a call to TPALT returns the static 
pressure P q and temperature T q at the altitude H in » The pressure is then 
expressed in lb per' sq ft and the density follows by the equation of state 
as 

RHO * P / [53. 32(32. 2) T ) (VI. 1.138) 

o o 

as expressed in slugs per cu ft. The Mach number to start the iteration 

is then computed as the maximum value, corresponding to the maximum 

operating airspeed V ; 

MO 

*MN " V MO /(29 ’ 0721 ^0 0023769 ) (VI. 1.139) 

The denominator of this expression is the speed of sound in kts, for 

any values of air density and temperature. Since • 0 , the first 

loop is entered in order to find the Mach number for minimum power 

required. The drag is initially ser equal to the weight, and the Mach 

number is reduced from the value X in steps of 0.1, until the value 

MN 

of drag returned by CTAER is greater than the value on the previous 
iteration. After this overshoot, the Mach number is increased by .01 
and on the next iteration, 

DRG > DRGp , 
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so control is transferred to statement 8, where X is the Mach 

MIN 

number at which drag is a minimum. 

The longer loop, between statements 4 and 10, begins by calling 

CTAER at the Mach number X. . The cruise drag returned is DRG, and 

MN 

set equal to the required thrust 


- DRG 

NRQ 


(VI. 1.140) 


before subroutine ENGINE is called, at the same values of Mach and alti- 
tude. The thrust per engine is returned, in terms of which the total 
thrust available is 


NAV 


EN * TH 
P PENG 


(VI. 1.141) 


where EN p is the number of engines, which is input in the canmon state- 
ment, UNIV. 

Now, if thrust available, F ... exceeds the thrust required, F„__ 

NAV NRQ 

on the first iteration (while J c ■ 0) it means that a velocity greater 

than the "maximum" value V is possible, and control passes to the 

MO 

write statement 20 before returning. Otherwise, the dimensionless 
error is 


err m “ (f nav ~ f nrq )/f nrq 


(VI. 1.142) 


and this difference, together with its previous value, permits use of a 
linear iteration procedure in subroutine ITRMHW, leading to a Mach 
number for which thrust available and thrust required are more nearly 
equal. If the error is less than .002, the computation is complete, 
and the true airspeed in kts is found at statement 10, 


V KTAS * - 5921 *HH * ‘’- 1 * ^ 


(VI.l f 143) 


following which the subroutine returns control to the calling program. 
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VI. 1.11 Subroutine RGBAL 


This subroutine is called only by MAIN, and Its purpose is to 

compare the total range or endurance computed in XRANGE with an input 

required mission range or endurance, R_ DD _. If R„ D _ <24, it assumes 

CRKy CRRQ 

an endurance requirement is being specified. 


RA ° I (1) • W * T CLHR <“ %RS < 24 -> 

If the error exceeds one percent, the aircraft gross weight is adjusted 
and the aircraft entirely resized before reflying the mission to 
generate an updated figure for range or endurance. This subroutine 
controls the program flow during this resizing iteration to determine 
the aircraft size required for the range or endurance requirement, 

R. . . „ . 

CRRQ 

A flow diagram for subroutine RGBAL is shown in' Figure Vl.l.t, which 
outlines the calculations corresponding to the following logic. When 
RGBAL is initially called, the range or endurance corresponding to 
the initial gross weight is known, as is the range or endurance 
requirement. Assuming the difference exceeds one percent, the gross 
weight is changed by a scale factor, FAC^. There are two options for 
resizing the aircraft, one holds wing loading fixed, the other fixes 
wing area, and is followed by reflying the mission and calculating the 
new range or endurance. This permits linear estimation of the gross 
weight needed to meet the required range or endurance condition. 
Recomputation leads to a new range or endurance and to a new estimate 
of the gross weight. A third iteration yields either less than l percent 
error or the printing of an error "failure" message before returning 
to MAIN. 

The subroutine begins with the initialization of the gross weight 
and range or endurance parameters, 


W GRI (1) " 


W_ 


(VI. 1.144) 


RAG I (1) 


r crnm + r clbnm 


(VI. 1.145) 


where 
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W G 

■ 

input gross weight, lb 

r crnm 

- 

cruise range, ran 

r clbnm 

- 

climb range, ran 

t crhx 

- 

time in cruise/hrs 

t clhr 

- 

time to climb, hr 8 


where the climb end crulee times on the right ere Input through 
COMMON/RANGX/, as computed in subroutine XRANGE. The scale factor 
FAC^ is input through C0MM0N/1NRBAL/, and the range error is expressed 
in dimensionless form as 


RC ERR “ (R CRNM + R CLBNM " r crrq )/r crrq 

or as 

rg err " <t clhr + t crhr ~ r crrq )/r crrq 


(VI. 1.146) 
(VI. 1.147) 


depending again upon whether or not Rcrrq Exceeds 24 in magnitude. 

If this figure is below .01 in magnitude, the subroutine returns, after 
printing the significant results. Otherwise, the flow depends on the 
iteration number which is initially 0. 

When ■ 0, the gross weight is increased or decreased by 

the scale factor, according to 


FAC U1 W c ’ 

rg err > 0 

(VI. 1.148) 

U C /FAC W1 • 

sg err <0 



and subsequently " 1 to indicate the first iteration is underway. 

For subsequent iterations, when ■ 1 or 2, gradient information 

is provided by 


SLOPE - ABS[(RG err - RG err1 )/(W c -W c1 )] 


(VI. 1.149) 
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and the up da tad gross weight is found as 


M c • M c - RC im ISt0Tt - 

after previously defining W C1 • Hg. 
The wing loading la redefined as 


(VX. 1.150) 


WGS " W G /S W 


(VI. 1.151) 


If 


^SWING “ and after some lines of output printing, control passes 


to sMtement 21. Subroutine SIZE is called twice, with NPC • 3 (to 
Initialize It) and NPC ■ 2 (to develop new values for wing and tail 
geometric parameters) . 


This is followed by three calls to subroutine FLAPS, corresponding 
to flap deflections in cruise, take-off and landing configurations. 

Lift and drag coefficients for the aircraft are returned for these flight 
conditions, and Mach number and altitude are defined as input to sub- 
routine CTAZR, which returns angle of attack, lift and drag in cruise. 


Subsequent calls to subroutine ENGSZ return values of engine sizing 
data, depending on flight condition and aircraft geometry. These are 

followed by calls to subroutines ENCWGT, WGHT, and PERFRM, which finally 

develops values for total range or endurance in the climb and cruise 

configuration. 

It is noted that at most three iterations are possible through the 
subroutine. If convergence occurs before the third iteration is complete, 
the variation of gross weight and corresponding total range or endurance 
is printed, together wi'th the required cruise range or endurance, R„ OBr .. 

Before returning, the subroutine calls CTAER, OUTPUT and AEROUT, 
to provide output of the final values of geometry, weight, and cruise, 
takeoff, and landing aerodynamics. 
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VZ.2 PERFORMANCE MODEL USER'S MANUAL 


As suggested by the descriptive material in Part X of Volume VI, 
the study of aircraft performance involves many subroutines and many 
input/output parameters. The performance subroutines given in Table 
VI. 1.1 are tabulated alphabetically in this section, and it will be noted 
that a relatively small number of input parameters (10 to 20) are 
required by each. 
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VI. 2.1 Input - Subroutine ACCEL 


VARIABLE 

DESCRIPTION 

EM 

■ initial Mach number 

ENP 

■ number of engines 

H 

■ altitude, ft 

SPEED 

■ desired velocity, kts, or (less than S) desired Mach number 

ST 

• initial time, hr 

SW 

■ wing area, sq ft 

M 

• 

■ initial gross weight, lb 


vi. 2. 2 Output - Subroutine ACCEL 

| VARIABLE 

DESCRIPTION 


■ Mach number 


■ range, nm 

ST 

■ time, hr 

W 

» gross weight, lb 

wr 

■ fuel consumed, lb 












2.3 Input - Subroutine ASPEED 


VARIABLE 

EM 

EMMO 

ENP 

H 

HIN 

KENGX 

KWRITE 

VMO 

WG 

HIN 


DESCRIPTION 

■ Mach number 

■ maximum operating Mach number 

■ number of engines 

■ altitude, ft 

• input altitude, ft 

■ indicates engine power setting 

0, do not write output 

1, write output 

» maximum operating velocity, KEAS 

* gross take-off weight, lb 

■ input weight, lb 



VARIABLE 

DRG 

XMN 


VI. 2.4 Output - Subroutine ASPEED 

DESCRIPTION 
= drag , lb 

= Mach number for cruise equilibrium 
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VI. 2. 5 Input - Subroutine CLIMB 


VARIABLE 

DESCRIPTION 

ALPHLO 

= angle of attack for zero lift, deg 

AR 

= aspect ratio 

DELH 

3 altitude increment in computing climb performance , ft 

DLMC4 

= quarter chord wing sweep, deg 

EM 

= Mach number 

EM 

CRU 

= cruise Mach number 

^MO 

= maximum operating Mach number 

EN P 

= number of engines 

EYE 

W 

= wing incidence angle, deg. 

H 

= altitude, ft 


= terminal altitude for climb maneuver, ft 

ENDCL 


*clm 

= 1, maximum ROC; 2, maximum operating speed; 3, specified EAS 

Sjrite 

= write option 

R 

= range, nm 

sw 

=> reference wing area, sq ft 

V 

CLMB 

=> climb velocity for maximum rate of climb, ft per sec 

V 

LI MX 

» limiting operating climb velocity 

V 

MO 

= maximum operating airspeed, kts 

W 

= initial weight, lb 
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Output - Subrout i ne CLIMB 


VARIABLE 

DESCRIPTION 

LM 

= Mach number 

II 

= altitude, ft 

R 

■ range , nm 

ST 

■ time, hr 

TA 

=» engine thrust, lb 

W 

= weiqht, lb 

WF 

=> fuel consumed, lb 
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VI. 2.7 Input - Subroutine DERIV 


VARIABLE 

DESCRIPTION 

ALPHLO 

■ zero angle of attack 

AR 

- aspect ratio 

CDFU 

• drag coefficient, take-off flaps, gear up 

CIALPH 

- lift curve slope, per deg 

CLMAX 

- maximum lift coefficient 

DLMC4 

• quarter chord sweep, deg 

DTFLP 

■ time increment to retract flaps, sec 

DTGR 

- time increment to retract gear, sec 

EM 

■ Mach number 

EYEW 

■ wing incidence, deg 

HOB 

m altitude above ground in wing spans 

IOUT 

to, all engines 
■ (l, one engine out 

12, accelerate-stop on runway 

MU 

- normal force friction coefficient 

QSW 

• product of wing area and dynamic pressure, lb 

RHO 

• air density, slug per cu ft 

TA 

• thrust, lb , 

TFLP 

■ time of beginning flap retraction, sec 

TG 

- time of beginning landing gear retraction, sec 

TGU 

■ time gear retraction completed, oec 

THEMAX 

» maximum fuselage pitch angle, deg 

VEND 

• velocity for maximum climb rate 

U 

■ weight, lb 

X1FMAX 

- maximum load factor during take-off 
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VI. 2 . 8 Output - Subroutine DERIV 

VARIABLE DESCRIPTION 

ALPHA = angle of attack, deg 

CD = drag coefficient during take-off 

CL => lift coefficient during take-off 

DCDF = drag coefficient increment due to flaps 

DCDG = drag coefficient increment due to gear 

DCLF => lift coefficient increment due to flaps 

UMAX =* altitude at end of take-off, ft 

IACCL = flag indicating limits on load factor or speed are reached 

THETAF = fuselage attitude angle during take-off, deg 


XLF 


load factor 






VI . 2 . 9 Input - Subroutine DLAND 


VARIABLE 

DESCRIPTION 

ALPHLO = 

! zero-lift angle of attack, deg 

ALTLND 

- altitude of runway above sea level, ft 

AR 

■ wing aspect ratio 

B 

■ wing span, ft 

CLALPH 

- lift curve slope, per rad 

DLMC4 

- quarter chord sweep angle, deg 

ENP 

- number of engines 

EYEW 

- wing Incidence angle, deg 

HAPP 

- approach altitude, ft 

HTC 

- wing height above ground, ft 

HER 

( 1, runway length constraint satisfied 
( 2, iterating on runway length constraint 


/ 1, first iteration in runway length computation 

IWGS 

- < 2, second iteration in runway length computation 


V 3, last Iteration in runway length computation 


/ 0, landing weight is gross weight 

IWLD 

■ \ 1, landing weight is weight at end of previous segment 


\ 2, landing weight is WLp CT times gross weight 

HUB 

- breaking coefficient of friction 

RSMX 

■ maximum rate of sink, ft per min 

SINKTD 

■ sink rate at touchdown, ft per sec 

SW 

- wing area, sq ft 

TDELAY 

- pilot's brake reaction time, sec 


(CONTINUED) 






VARIABLE 


DESCRIPTION 


TDELLD 

o 

= incremental temperature above standard temperature, F 

TIDLE 

= idle thrust, lb 

TROT ID 

= ratio of reverse thrust to idle thrust 

VRATT 

= ratio of approach speed to stall speed 

W 

= aircraft weight, lb 

WG 

= gross weight, lb 

WGS 

=* wing loading, lb per sq ft. 

WLPCT 

= landing weight divided by gross weight 

XLDGRQ 

= landing length requirement, ft 

XLFMX 

= maximum load factor during flare maneuver 


VI. 2. 10 Output - Subroutine DLAND 


VARIABLE 

DESCRIPTION 

TIDLE 

= total idle thrust, lb 

WGS 

= wing loading, lb per sq ft 

DLT 

= total landing distance from obstacle height, ft 

DLG 

= distance from obstacle to start of flare, ft 

DI/TR 

= flare distance, ft 

DDE LAY 

= distance from brake application, ft 

DLGL 

= braked ground roll distance, ft 

DFAR 

= FAR 25/121 factored landing distance, ft 

RS 

= rate of sink during glide portion, ft per min 

VAPEAS 

= approach speed, KEAS 

VSTEAS 

= landing stall speed, KEAS 

VTD 

= touchdown speed, KEAS 

DOG 

= average decceleration during landing ground roll, g's 

THETA 

= approach flight path angle, deg 

Cl MX 

= landing maximum lift coefficient 
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VI. 2. 11 Input - Subroutine PERFRM 


VARIABLE 

DESCRIPTION 

CRALT 

= cruise altitude, ft 

EMTURN 

* Mach number in turn 

HPORT 

= altitude of airport, ft 

HTURN 

= altitude of turn, ft 

I DC 

* control parameter 

IFLY 

= input parameter to control mission computation 

ISEG 

= 1 to 8, mission segment indicator 

MUB 

= braking coefficient of friction 

TDELTO 

= temperature increment above standard for mission 

TDELTX 

= temperature increment above standard for engine 


sizing 

VLIMX 

= maximum take-off velocity, ft per sec 

WG 

= aircraft gross weight, lb 

WTMISN 

= start of mission weight, lb 
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VI. 2. 12 


Input - Subroutine PGBAL 


VARIABLE 

DESCRIPTION 

DFLPLD 

■ landing flap deflection, deg 

DFLPTO 

- take-off flap deflection, deg 

EMCRU 

■ cruise Mach number 

FACW1 

■ weight modification factor 

HNCRU 

• cruise altitude, ft 

IEGWGT 

i 0, propulsion weight and nacelle geomtry fixed input 
' 1, propulsion weight and nacelle geomtry computed in ENGWGT 

ISWING 

m t 0, hold wing loading fixed 


* 1, hold wing area fixed 

NT YE 

■ input flag l to 13, describing engine type 

NTYP 

■ input flag 1 to 16, describing propeller type 

RCLBNM 

• climb range, nautical miles 

RCRNM 

■ cruise range, nautical miles 

RCRRQ 

- required total mission range, nautical miles or endurance, hrs 

SW 

- wing area, sq ft 

TCLHR 

■ climb time, hrs 

TCRHR 

- cruise time, hrs 

WG 

■ initial gross weight, lb 

VI. 2. 13 OUtput - Subroutine RGBAL 

VARIABLE 

DESCRIPTION 

WG 

■ final gross weight required to meet range condition, ib 

WGS 

- final wing loading, lb per sq ft 
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Vi. 2. 14 Input - Subroutine TAKOPF 


VARIABLE 

DESCRIPTION 

ALPHLO 

■ zero lift angle of attack, deg 

AR 

■ aspect ratio 

B 

■ wing span, ft 

CLALPH 

■ lift curve slope, per rad 

CLMAX 

■ maximum lift coefficient at stall 

DELTVR 

* time required to rotate fuselage in take-off, sec. 

DLMC4 

■ quarter chord sweep angle, deg 

DVR 

- increment of rotation speed above decision speed, kts 

DV1 

■ increment of decision speed above stall speed, kts 

EYEW 

■ wing incidence angle, deg 

CRCD 

■ drag coefficient increment of landing gear 

HFLP 

- altitude at which flaps are retracted, ft 

HTMAX 

« maximum take-off altitude, ft 

MUB 

• Braking coefficient of friction 

NFA1L 

| 0, one engine is out 
( 1, all engines operative 

UM 

- rolling coefficient of friction 

VCLMB 

■ velocity for best rate of climb, kts 

VRAT 

= margin of take-off speed above stall speed (greater 
than one) 

W 

= current aircraft weight, lb 

WF 

= current fuel weight, lb 

WGS 

= wing loading, lb per sq ft 

XLF 

= load factor during takeoff 

XTORQ 

= required field length, ft 
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VI. 2. 15 Output - Subroutine TAKOFP 

VARIABLE 

DESCRIPTION 

ALPHA 

• angle of attack, deg 

CD 

• drag coefficient 

CL 

■ lift coefficient 

DTFLP 

- time required to retract flaps, sec 

DTCR 

“ time required to retract landing gear, aec. 

H 

• altitude during take-off, ft 

HMAX 

• cutoff maximum altitude, ft 

QSW 

* product of dynamic pressure and vlng area, lb 

RELI 

■ Reynolds number per ft of length 

RHO 

- air density, slugs per cu ft 

ST 

■ time end of take-off 

TA 

• total engine thrust, lb 

TFLP 

■ time at start of flap retraction, sec 

TFLPU 

~ time at end of flap retraction, dec 

TG 

■ ti<*e of start of gear retraction, sec 

TGU 

* time at end of grear retraction, sec 

VEND 

» best rate of climb speed, kts 

VLLMX 

■ constrained rate of climb speed, kts 

W 

- aircraft weight, lb 

WF 

- weight of fuel consumed, lb 

XTO 

■ ell-engine runway distance needed to clear 35 ft 


altitude, ft 

DSTOl' 

= accelerate-stop distance, ft 

S35 

= engine-out distance to 35 ft altitude, ft 

S35ALL 

= all-engine distance to 35 ft altitude, ft 

SFAR25 

= FAR 25 factored takeoff distance, ft 

ssn 

- all-engine distance to 50 ft altitude, ft 
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VI. 2. 16 Input - Subroutine TAXI 


VARIABLE 

DESCRIPTION 

DELTT 

■ time Interval, hr 

ENP 

- number o£ engines 

H 

■ altitude, ft 

KWRITE 

| 0, do not write output 

m / 

( 1, write output 

ST 

- initial time, hr 

TA 

» engine thrust, lb 

W 

- initial aircraft weight, lb 

WT 

- initial burned fuel weight, lb 


VI. 2 . 17 Output - Subroutine TAXI 


VARIABLE 

DESCRIPTION 

ST 

» final time, hr 

W 

- final aircraft weight, lb 

V? 

■ final burned fuel weight, lb 


62 










VI. 2. 18 Input - Subroutine TURN 


VARIABLE 

DESCRIPTION 

CLTLMT 

■ maximum lift coefficient in turn 

EM 

- Mach number in turn 

EN P 

- number of engines 

HTURN 

■ altitude of steady turn, ft 

I DC 

» performance or engine sizing control parameter 

SW 

« wing area, sq ft 

W 

■ aircraft weight, lb 

XLFTRN 

* desired load factor in turn 


VI. 2.19 Output - Subroutine TURN 


VARIABLE 

DESCRIPTION 

TRTKN 

XLFTRN 

thrust required in steady turn, lb 
load factor in turn 
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VI. 2.20 Input - Subroutine XRANGE 


VARIABLE 

DESCRIPTION 

CRMACH 

- cruise Mach number 

DMACH 

* wing compressibility Mach number 

EM 

* Mach number at end of climb 

EMMO 

-* maximum operating Mach number 

ENP 

■ number of engines 

EYEW 

» wing incidence, deg 

FRES!' 

* reser/e fuel fraction (1.0 equals 45 min> 


(greater 10. equal*, lbs 

H 

=■ altitude, ft fuel) 

HNCRU 

= design cruise altitude, ft 

HPAVLB 

= horsepower available 

KODE 

= engine operating condition indicator 

SUP 

= shaft horsepower 

ST 

= time to begin cruise, sec 

VMO 

= maximum operating air speed, kts 

W 

= weight, lb 

WF 

= fuel weight, lb 

WFAMAX 

= maximum weight of available fuel, lb 

WFAMIN 

= minimum weight of available fuel, lb 

WG 

= gross weight', lb 

WTMISN 

= weight at start of mission, lb 

ZLQD 

= lift to drag ratio 
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VI. 2.21 Output - Subroutine X RANGE 


VARIABLE 

DESCRIPTION 

EM 

= Mach number 

EMCOST 

* Mach number used in cost calculations 

EMSAVE 

^ Mach number used in RGBAL calculations 

FPMN 

= Mach number used in fixed pitch propeller iteration 

H 

= altitude, ft 

OUT (6, 16 

) flight conditions at start and end of cruise, for three 


Mach numbers 

PCRPM 

* percent of rated rpm 

PO 

= static pressure, lb per sq ft 

R 

= total range, nm 

RCLBNM 

= climb range, nm 

RCRNM 

= cruise range, nm 

ST 

= total time, hr 

TCLHR 

= climb time, hr 

TCRHR 

- cruise time, hr 

W 

= weight, lb 

WF 

= fuel weight, lb 

WECRU 

= cruise fuel weight, lb 

WFRES 

= reserve fuel weight, lb 
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VI. 3 PERFORMANCE MODEL PROGRAMMER'S MANUAL 


The flow charts in the following pages to,; low the order shown in 
Table VI. 1.1, and it will be noted that the nominal mission sequence 
of Figure VI. 1.1 makes iterative loops largely unnecessary. That is, 
each mission segment requires calling a different subroutine concerned 
with the segment, and the output of each segment acts as the input to 
the next. The individual segments are air j straghtforward , although 
the flight condition must satisfy variou constraints, as tested by a 
sequence of IF-tests. 
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VI . 3 . 1 Subroutine PERFRM, Mission Computation Control Routine 


Subroutine PERFRM controls the mission computations in GASP, 
of mission segments coverings 

1. Taxi 

2. Take-off 

3. Acceleration 

4. Climb 

5. Range 

6. Landing 

7. Turns 


A series 


are controlled by PERFRM. The routine PERFRM is called from the main 
program. 
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PERFRM 


2 






















VI. 3.2 Subroutine ACCEL, Acceleration Segment 


This routine computes level flight accelerations under the control 
of PERFRM. Several subroutines are called by ACCEL during these calcula- 
tions including: 

TPALT — Atmospheric Properties 

DRAG — Drag Calculations 

ENGINE — Propulsion Characteristics 

Major outputs computed by this routine are weight, range and time changes 
during the acceleration. 
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ACCEL 



Figure VI. 3. 2 Subroutine ACCEL, Detailed Flow Chart 
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ACCEL 


2 


T=FN*ENP 
WFUEL=WFUEL*ENP 
TAS=EM*SA/1 . 69 
EAS=TAS*SQRT (RHO/ . 0023769) 


c 


ic 


X 


KWRITE 


>0 


WRITE(6)ST,R, 

WF,W,H,TAS,EAS, 

EM.EMD.T.WFUEL 







VI. 3.3 Subroutine CLIMB, Climb Calculations 


This routine computes planar climb segment characteristics. Atmospheric 
properties along the climb path are determined from TP ALT. Required lift 
coefficient is found from CLIFT and drag coefficient from DRAG. Propulsion 
system characteristics are obtained from ENGINE. 

Climb paths may be flown at maximum rate-of-climb, at maximum allowable 
airspeed, at a specified equivalent air speed. Three functions PART, 

BOZA, and ADEN are defined and used in subroutine CLIMB. 
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CLIMB 


FUNCTION DEFINITIONS 


PART(XXX , X 1 1 1 , X222 , X? , . : 1 1) = ( ((XXX-X11 1) * (XXX-X222) )/ 
((X333-X111)* (X333-X222 ) ) ) # Y1 1 1 
B02A(V1 1,V22,V33,U1,U2,U3)=(V 11**2* (U2--U3) -V22**2* (U1 -U3) 

♦V33**2*(Ul-U2))/(2.*(-Vll*(U2-U3)*V22*(Ul-U3)-V33*(Ul-U2) 
ADEN=W1 , W2 ,W3, 21 , Z2 , Z3) =Z3* (W2-W1 ) *22* (W1 -W3)+Z1* (W3-W2) 


RTOl)=57. 2957795 

1NDIOO 

deltaino. 

ICLMP*1CLM 
ALPIIA= 12. 


CALI. CI.I IT( 1, ALPHA, CLL,ALPHO,CLALPH, HOB, AR.DLMC4, EM) 


WRITE (6) 
MEN DC L 


WRITE (6) 


WRITE (6) 
ilENDCL 


CALL TPALT (II , GALT , PO , FKALT , TO , G , XKV) 


P0=P0* 144 . 

SA=49. 1 *SQRT (TO) / 1 . 689 
RIK)=P0/53. 32/TO/G 
TAMB=T0 
PAM B-PO 

VL=1. 1 0*SQRT (W/ (SW* 1 . 42636*RHO*CLL) ) 


LC 1-0 
l.C2=0 
l.C3=0 
KKK=0 

1CIA1= ICIMP 

VMAX VMO/SQRT(R!IO/ . 0023769) 


Fiyuie VI. 3. 3 Subroutine CLIMB, Detailed Flow Chart 
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CLIMB 


3 
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CLIMB 


5 


RCTA=PART (V2,V22,V33,V11 , RCT1)+PART (V2,V11 ,V33,V22,RCT2) 
♦PART (V2,V11,V22,V33, RCT3) 

RCPA=PART (V2,V22,V33,V11, RCP1) ♦PART (V2 , Vll , V33, V22,RCT2) 
♦PART(V2,V11 , V22,V33,RCP3) 
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CLIMB 


6 















CLIMB 


7 










VI. 3. 4 Subroutine DLAND, Landing Calculations 


Subroutine DLAND computes landing segment characteristics from an 
analysis of terminal glide, flare, touchdown and braking roll measures. 
Atmospheric properties for the landing field are computed from TPALT . 
Routine CLIFT and DRAG provide landing aerodynamic characteristics, 
propulsion system characteristics are obtained through the routine ENGINE. 


01 


VRAT=VRATT 

TIDl.E=0. 


DLAND 


>* 99999 . 



WKITT.(6')Wf,SM\X 


Figure VI. 3. 4 Subroutine DLAND, Detailed Flew Chart 
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I 














3 


DLGL=HAPP/TAN (THETA) 

DLTR= ( (RZ*THHTA) 12 . )* ((1 . O-CAMTD/THETA) 


VTP-(VRATM.‘)*.5*VST 

l)ni.l.AV--V"rn*TI)I:LAY 

HOBIITR/B 

! Ml - VTI)/SQRT( 1 .4*0*53. 32*TAMB) 


t'M.I. ai IT (1 , 1.YI.W, GI.RL, ALl’IILO ,CLAPH,HOB, AR.DLMC4 , EML) 
GM.l DRAG ( 7 , EMI, , HOB . CI.RI, , EMU , CURL) 


PLRL-CDRL- (MlJB*CllRi,l 

ARAT- PI.R1./ (Cl. MX* (VST/VTP) **2 ] 

TOW!.- -T TUEE* TRO P! !>/ WI. 


TOWL= -100. *TROT 1 D/ WL 


TROTH) 


l ; R\T=Mlli! FOWL 

Al.N- A LOG ( ! : RAT / ( L RAT +A RAT) ) 

PLG-- 1 3. 0287*WOS*ALN'/ ( DRAT*DLRL) 
DOG - VTP* * 2. / ( Dl,G*2 . *G) 

DI. r=DlA>Dl.TR + PDI'l AY-t-Dl Gi. 


XLIH.RO 


= 01)111)1). 


R L I URN 









DLAND 


4 



I 


RETURN 



ii.i:r=i 

TH1TA-TI1HTA*57 . 29 
V(il,-V('.I,/ 1 .689 
VAIM:AS=V(iF,*SRDRAT 
VSTi:AS*VSTEAS/ 1.689 
m ; \R-m.T/.6Q 
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VI. 3. 5 Subroutine TAXOFF , Take-off Calculations 


Subroutine TAKOFF computes take-off performance through ground roll, 
rotation, lift-off, gear retraction, initial climb, and flap retraction. 
The take-off trajectory is integrated and derivatives are computed by 
subroutine DERIV. Atmospheric properties are computed by TPALT. 
Aerodynamics are obtained from CLIFT and DRAG. Propulsion characteristics 
are obtained from ENGINE. Integration is performed by the utility routine 
1NTS . 


86 



TAKOFF 


I YYY (XXX.X1,X2,Y1.Y2)~Y1»(XXX-X1)*(Y2-Y1)/(X2-X1) I 


UMAX=1ITMAX*H 

WSTO-W 

WFSTO=W 

ioirr=o 


FUNCTION 

DEFINITION 


TG=999. 

TFLP=999. 

MU=UM 
M= 0 . 

I ACC 1=0. 

I M= 0 . 

IR=0. 

I PAS=0. 

Xl,0ST=l . 

1 1* L= 0 . 

IG=0. 

IK'IKt=0. 

!H'!)F=0. 

ih:i.f=o. 

MC,R = 50. 
mr>R=7. 

III.!- P-400 . 

m FI, P--PLFPTO/3. 33333 
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TAKOFF 


5 


TIM1*YYY (VI ,PEAS,EAS,PT(2) ,T(2)) 
VW= YYY (V 1 , PEAS , EAS , PVKTS , VKTS ) 
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TAKOFF 


7 


T.4MB-T0+TDELT 

P0*P0*144. 

SA*49. 1*SQRT(TAMB) 
RHO*PO/53. 32/TAMB/G 
RELI*T(4)/XKV 
PAMB^PO 

HOB a (T(7)4HTG-HAPT)/B 
EM-T (4)/SA 
VKTS=T(4)*. 592087 
EAS=VKTS*SQRT (RHO/ . 0023769 
QSW».5*RHO*SW*T(4)**2 
KENG=5 


t 



♦ 


l THETA 

t F»Tf5)«57.. 


> 

, I 


| T(1)=T1SAV 



1 CALL I NTS (T, 4, 2 


i 




□EEJ 


c 


KWRITE 




k. 



(225 
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TAKOFF 


8 
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TAKOFF 


9 


HHH-50.+HAPT 

S50-YYY(HHH,PT(7) ,T(7) ,PT(6) ,T(6) 
I PAS* 2 
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TAKOFF 


10 
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VI. 3.6 Subroutine DERIV, Tine Derivative Calculation 

Subroutine DERIV computes the instantaneous flight path and velocity 
derivatives required by subroutine TAKOFF. Routines CLIFT and DRAG 
provide aerodynamic character! sties. Propulsion system characteristics 
are computed from ENGINE. 
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DERIV 
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DERIV 


2 
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DERIV 


3 










c 


DERIV 


T (9) - (32 . lb/W*T (4) ) ) * (CL*QSW*TA*SIN (ALPHAR) -W**COS (T (S) ) ) 
T ( 10) «T (4) *C0S (T (5) ) 

T(11)«T(4)*SIN(T(S)) 


RETURN 


< 


( 
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VI. 3.7 S ubroutine TAXI. Taxi CaIcul r »ton « 


Subrouting TAXI conputgs taxi fuel requirgsjgnts . Subrouting TPAW 
guppligg atmospheric conditions at thg runway altitude. Subroutine ENGINE 
supplies propulsion system char ac tar istics during ground idle. 



RETURN 


Figure VX.3.7 Subroutine TAXI, Detailed Flow Chart 
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VI. 3.8 Subroutine TURN. Turning Performanc e 


Subroutine TURN computes steady state turn maneuvers. The routine 
TP alt provides atmospheric properties during the turning maneuver. Drag 
coefficient is obtained from DRAG and propulsion system characteristics 
from ENGINE. 



RETURN 


WRITE (6) 
TA.TRTRN 



Figure VI. 3.8 Subroutine TURN, Detailed Flow Chart 
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VI. 3. 9 Subroutine X RANGE, Cruise Range Calculation 

Subroutine X RANGE computes cruise performance in GASP. Calculations 
are based on the Brequet equation. Atmospheric properties are computed 
by TP ALT. Maximum speed is obtained from ASPEED. Cruise aerodynamics 
are computed by CTAER and propulsion characteristics from ENGINE. 
Accelerating flight paths are computed by ACCEL. 
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XRANGE 


KWRITI-KWRITE 

ST1-ST 

NOFUEL* 0 

R1»R 

WF1=WF 

W1*W 

DELWGT-WG- (W1+WF1) 

HN1*H 

INDFP=0 


, \ — — - 

|CALLTPALTOj^ 


PO=PO*144. 

ASON=49. 1*SQRT(T0) 

EMX=EM 

SQRDR=SQRT (P0*S18 . 67/(2116. 22*T0) ) 

1 


c 


KODE 


ISEGX=0 

KENG=3 


*8 


c 


NTYE 


I 


» Ikeng.5 I 




SPECR=0. 

DELV=-10. 

ITOM=0 

VELK=EMNP*29.07*SQRT(T0) 


| EM=VELK/29. 07/SQRT(T0) | 


3 


EM 


>EMNP 


1 

1 

1 SPECRP=SPECR 1 

1 

1 


C 


KWRITI 


> 


I CALL CTAeH 


-|kwrite*o. 1 


|Tnrq=drg/enpT 

I 

Figure VI. 3. 9 Subroutine XRANGE, Detailed Flow Chart 
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XRANGE 


2 



i 


— ►DO 105 J*l, 16 


| OUT (I ,J)=0. | 
L< * i 105 continue! 
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XRANGE 


4 


KSOKESa 



KWRITE=KWRITI 

ENRQ=DRG/ENP 



VKTS1=VK*.S921 

VKES1=VKTS*SQRDR 

CL1=CL 

EM1=SMN 

EMD1=DMACH 

ALPlil=AD 
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Ill 









XRANGE 


6 


WRITE (6) 
SMN.JC 



FF2-FF*ENP 

INDFP-0 

EM2*SMN 

VX-EM2*AS0N 

VKTS2«VX*.5921 

VKES2-VKTS2*SQRDR 

BREG2«VX*W283600. /FF2/6076. 1 

BREGA- (BREGl+BREG2)/2. 

RCRNM- BREGA* ALOG (1 . / (1 . -WFCRU/W1) ) 

R2-R1+RCRNM 

EMD2-DMACH 

ALPH2-AD 

THET2-AD-EYEW 

CL2=CL 

ZLOD2-ZLQD 

VAVG= (VKTSl*VKTS2)/2. 
ST2-ST1+RCRNM/VAVG 
KK=2*K-1 
OUT (KK, 1)=ST1 
OUT(KK,2)=Rl 
OUT (KK, 3)=WF1 
OUT(KK,4)*Wl 
OUT(KK,5)=HNl 
OUT(KK.6)=VKTSl 
OUT (KK, 7)=VKES1 
OUT (KK, 8) -EMI 
0UT(KK,9)-EMD1 
OUT ( KK , 1 0) -ALPH1 
OUT(KK, 11J-THET1 
OUT(KK, 12)=CL1 
OUT(KK, 13)=ZL0D1 
OUT (KK, 14)-FF1 
OUT(KK, 15)=BREG1 
OUT(KK, 16)=JDUT(KK,6)/OUT(KK, 14) 
OUT(KK«l,l)=ST2 
OUT (KK+1 ,2)-R2 
OUT(KK+l,3)=WF2 
OUT (KK+1 , 4)=W2 
OUT (KK+1 ,5)=HN1 
OUT (KK+ 1 , 6) =VKTS2 
OUT (KK+1 , 7)=VKES2 
OUT(KK+l,8)=EM2 
OUT (KK+1 , 9)=EMD2 
OUT (KK+1 , 10)=ALPH2 
OUT ( KK+ 1 , 1 1 ) -THET2 
OUT(KK+l, 12)=CL2 
OUT (KK+1 , 13) =ZLOD2 
OUT(KK+l, 14)=FF2 
OUT (KK+1, 15) = BREG2 
OUT(KK+l , 16)=OUT(KK+l,6)/ 
0UT(KK+1,14) 

RESF (K)-WFRES 












XRANGE 


7 





^ 90 ^- 


WFAMAX 


J 




WFAVLB 


j WFAVLB-WFAMAX j - 


-0 


c 


WFADES 


■ 


WFAVLB 


I WFAVLB*WFADES \ - 

^ I 

j WRITE(6)j 

| NOFUEL" 1 1 

♦ 


-0 


c 


WFAMAX 


> 


>WFAVLB 



250] 


KWRITE*KWRITI 
KK=2* (ICUS+1) -1 
RC L BNM=OUT ( KK , 2 ) 

RCRNM=OUT ( KK* 1,2) -OUT ( KK , 2 ) 

TCRHR= (OUT (KK*1 , 2) -RCLBNM) /OUT (KK,6) 
TC LHR= (OUT ( KK , 1 ) 

EMSAVE=OUT(l,8) 


c 


J. 


KODE 


> 


=8 


" |EMSAVE=OUT (KK,8) j 


HSAVE=OUT ( K K , 5 ) 

H=HSAVE 

WFCRU=OUT(KK+ 1,3) -OUT (KK , 3) 

WFRES=RESF(ICRUS*1) 

R=OUT(KK*l,2) 

ST=OUT(KK+r, 1) 
WF=OUT(KK*l,3) 



RETURN 


U3 




I 

! ? f 

f*; 

VI. 3. 10 Subroutine RGBAL Range Balancing Calculation 

I Subroutine RGBAL balances the range computed by XRANGE against a 

desired Input range value by adjusting the take-off gross weight. The 
following subroutines are required in this resizing procedure: 

SIZE— Aircraft Geometry 

( FLAPS , CTAER— -Aerodynamics 

engine, engwt — E ngine Performance Size 
WGHT — Weights 
PERFRM — Performance 

( XRANGE — Cruise Range 
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Figure VI. 3. 10 Subroutine RGBAL, Detailed Flow Chut 
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RGBAL 


1 CALL WGHT I 

|| 

■ 

ICOND-O 

ISEG*0 

1 


I CALL PERFRM(ISEG.ICOND?IFLY!ni 

<f 


WRITE(6)RGERR, 
RGERR1,WG,WG1, SLOPE 



EM=EM 

H=UNC 

Bl 


■A 

fcALL 
I CALL 
1 CALL 

CTAER | 
OUTPUT I 
AEROUTl 


RETURN 


S 

\ 
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